Recoil spectrometers for heavy-ion 
identification and secondary-beam production: 
pushing the low-energy limit 

L. Audouin L. Tassan-Got^, P. Armbmster ^, 
K.-H. Schmidt ^ C. Stephan^, J. Taieb^''^''^ 

^FZK, Hermann-von-Helmotz Platz 1, 76344 Ettlingen-Leopoldshafen, Germany 
^IPN Orsay, CNRS-IN2P3, Campus Paris XI hat. 102, 91406 Orsay, France 
^GSI, Planckstrasse 1, 64291 Darmstadt, Germany 
•^DEN, CEA Saclay, 91191 Gif-sur-Yvette, France 



Abstract 

The feasibility of low-energy fragmentation experiments using a magnetic spectrom- 
eter is discussed. The main challenge is the multiplicity of the ionic charge states, 
which can hamper the identification in both Z and A of the fragments. Three topics 
are covered. First, a specific set-up for ionization chambers, based on a very large 
gas thickness, is presented. Its satisfactory performances are discussed in light of the 
observations during a 500A MeV Pb-|-p experiment performed at the FRS (GSI). 
As a second topic, the possibility to use a thick layer of matter (a degrader) as 
a passive measurement device to identify the nuclear charge and the ionic charge 
state of fragments is discussed. This method, successfully used for Z identification 
in experiments such as Pb-|-p at lA GeV, fails to measure the charge states at 
500A MeV for the same system. It is shown that surface defects of the degrader 
are probably responsible for this failure. The third topic is the description of new 
analysis techniques developed in order to account for and subtract the contribution 
of polluting charge states in the spectrometer, thus making possible a clean estima- 
tion of the production cross sections of all fragments. The combination of those new 
experimental and analysis techniques made the 500A MeV spallation experiment a 
success. 
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Introduction 



An in-flight fragment separation apparatus can be used in two ways: 

• as a recoil spectrometer, to measure production cross sections or kinematics 
properties of the fragments; 

• as a separation device, aiming at producing high-purity secondary beams of 
nuchdes. 

The Fragment Recoil Separator (FRS) [1] of GSI has been very successfully 
used for both purposes for more than ten years. By now, several projects 
across the world are dedicated to the construction of more powerful devices, 
in terms of larger fragment acceptance (spectrometer for fission products), 
and/or increased selectivity and high- intensity secondary-beam production. 

The upper energy range covered by such devices seems to be mostly a techno- 
logical (and financial) question, as increasing the energy requires more pow- 
erful magnets and a longer flight path, not even talking about the accelerator 
itself. But the lower range is determined by physics: it is assumed that only a 
very low proportion of non-fuUy stripped fragments can be tolerated in order 
to obtain a clean separation of fragments and to measure production cross sec- 
tions. Therefore, a whole energy range is somehow forbidden to fragmentation 
experiments. 

This low-energy limit has been faced by a recent experiment conducted at the 
GSI spectrometer FRS. The aim of this experiment was to use the inverse- 
kinematics method to measure production cross sections from residues of the 
spallation of lead by protons with an incident energy of 500A MeV. The results 
of this experiment will be the subject of a separate, forthcoming publication. 
This paper is dedicated to the description of the experimental and analysis 
techniques which were developed for this experiment and to discuss whether 
or not they proved to be successful. Part of these techniques may be used more 
generally than for the FRS setup from GSI. 

First we will briefly remind the main characteristics of the FRS (section 1). 
The Fragment Recoil Separator and its set of detectors have been designed 
using the hypothesis that most, if not all the fragment one wishes to study are 
fully stripped [1] . This is of course less and less true with decreasing energy. In 
section 2 we will discuss the variations of the charge-state probabilities with 
energy and their consequences on the selectivity of the spectrometer, as well 
as the possible choices of so-called stripping materials. 

With a lead beam energy as low as 500A MeV, no stripper is efficient enough 
so that the fraction of non-fuUy stripped fragments would be negligible. The 
response of the ionization chambers (used for the determination of the atomic 
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number of the fragments) in this situation wiU be discussed in section 3, as 
well as a new setup, developed for this experiment, based on an increased gas 
thickness. 

In section 4 we will describe a method to use a thick matter of layer as a 
passive measurement device in order to estimate the atomic number and the 
charge state of the fragments. We will also discuss, how unexpected problems 
prevented the measurement of the ionic charge state of ions (and thus of the 
mass) in the above-mentioned experiment. The final sections will be dedicated 
to new analysis techniques aiming at obtaining production cross sections, even 
in the case of an incomplete event-by-event mass identification. Section 5 will 
present a method to determine the most-likely mass of each fragment and to 
reduce the number of possible contaminating charge states. Then, in section 6, 
we will explain how to account for the large fraction of non fully-stripped ions 
and how to obtain isotopic production rates with low uncertainties despite the 
ambiguities in the mass determination. 

We will then conclude and try to answer whether or not the widely admitted 
low-energy limit has been pushed downwards. 



1 Main characteristics of the Fragment Recoil Separator 

The FRS, which is schematically presented in figure 1, has been designed to 
separate and identify nuclei produced by fragmentation of relativistic nuclei, 
the reference beam being lA GeV U [1]. Basically it is designed as a two- 
stage device: fragments produced in a target located at the entrance of the 
spectrometer (5*0) are selected according to their magnetic rigidity, reach the 
dispersive focal plane {S2), are then selected again in a second magnetic section 
and reach the final focal plane (5*4). The first magnetic section is intended to 
measure the magnetic rigidity of the fragments right after the target as well 
as to remove the primary beam in order to protect the detectors located at 
5*2 and S4. 

The behavior of a charged particle of mass m and charge Q moving in a 
magnetic field B is given by the well-known relation : 



Here p is the curvature radius of the particle trajectory in the magnetic field, 
c is the speed of light and /3 and 7 are the usual Lorentz relativistic factors. 



Bp^ 
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Fig. 1. Schematic view of the Fragment Recoil Separator. Each magnetic section 
between the focal planes consists of two dipoles (large blocks) plus several quadrupoles 
and sextupoles (arrows). The horizontal positions at S2 and S4 are usually measured 
using plastic scintillators, which are also used for time- of- flight measurement. 



For a nucleus with A nucleons, equation 1 can be reduced to: 
e.c q 



(2) 



u being the elementary nuclear mass, e the elementary electric charge, and 
q the number of ionic charges of the considered ion. In the rest of this pa- 
per, for simplicity, we will refer to the magnetic rigidity of a fragment as a 
dimensionless term, A/q.P'j. 

In experiments at the FRS, the A/q ratio of each nucleus is deduced from 
the measurement of its time of flight and magnetic rigidity in the second 
magnetic section, using plastic scintillators [2] . Z is deduced from energy loss 
in ionization chambers set at the exit of the spectrometer. 

One of the main hypotheses used for the design of the FRS was that most 
of the fragments should be fully stripped during their flight in the magnetic 
sections. In this case one has q = Z and getting the mass of each fragment is 
straightforward. 

As will be outlined in section 2.3, the insertion of a thick layer of matter at 
5*2 is a way to drastically reduce the number of fragments transmitted down 
to the exit of the spectrometer. Detailed calculations demonstrating this effect 
can be found in [3]. Let us point out that the possibility to choose the exact 
profile of the degrader allows one to use the FRS in particular magnetic optics 
modes: achromatic or mono-energetic. 

• achromatic: the position of the fragments at 6*4 does not depend on their 
velocity at the entrance of the spectrometer. This is particulary helpful to 
achieve a better separation of the fragments and was the mode selected for 
the 500 A Me V experiment which is the reference for this paper. 

• mono-energetic: the velocity of the fragments at S4 does not depend on their 
velocity at the entrance of the spectrometer. This mode is mainly used for 
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implantation experiments. 

In the following sections, we will often take the ^°^Pb+p at 500A MeV experi- 
ment as a reference for our discussions. In this experiment, the target consisted 
of 87 mg.cm"^ liquid hydrogen enclosed in 4 Ti foils, each 9 mg.cm^^ thick 
(see [4] for a complete description). After this target, a stripper foil of Nb 
with a thickness of 60 mg.cm^^ was set. At 5*2 the plastic scintillator had a 
thickness of 3 mm and was followed by a profiled degrader with a mean thick- 
ness of 1740 mg.cm"^ of pure Al. Considering the energy of the fragments 
after the degrader (roughly 300A MeV) , a Ti foil would have been suitable for 
optimal stripping, but Al was a good choice, nevertheless (see chapter 2.2 for 
a discussion about stripper efficiency). 



2 Physics of atomic charge states and consequences 

From equation 2 it is obvious that a magnetic separation device is sensitive to 
the ionic charge of the nuclei. This is also the case of energy losses (which are 
the key process for many detectors), as the elemental energy loss in matter is 
governed by electromagnetic interactions and therefore depends, not on Z'^ as 
one often thinks in the first place, but on q^ff, the square of the effective ionic 
charge of the nuclei. In the case of a gas medium, this effective charge state is 
almost equal to the charge state as the time interval between two interactions 
of the fragment with gas atoms is much larger than the decay time of atomic 
excited states of the fragment. 

2. 1 The physics of atomic charge states: a brief overview 

The evolution of the ionic charge state of an ion in matter depends on two 
competing processes: electron stripping and electron capture. In the relativistic 
or quasi-relativistic regime, the stripping probability cr^ is nearly independent 
of the energy of the ion. On the other hand, the probability of electron capture 
leading to the presence of the electron on shell X, a^^\ strongly increases with 
decreasing energy (for detailed discussion of the role of the various atomic 
shells and related phenomena such as the de-excitation time of the capturing 
ion, see for example [5]). One can distinguish two extreme situations: 

• a^^^ « as- any electron carried by the ion is removed after a very short dis- 
tance in the medium. Similarly, any electron capture is quickly followed by 
the stripping of this electron. Therefore the ion is mostly fully stripped. This 
situation corresponds to the high-energy limit (roughly, beyond lA GeV for 
heavy ions such as lead). 
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• cr^^) >> cr^: any electron stripped for the X shell is quickly replaced. This 
is most likely to happen if the velocity of the ion is close to the one the 
electron would have on its shell. In this situation the ion has its X shell 
filled most of the time. For example, considering heavy ions such as lead, 
this situation appears in the case of the K shell around lOOA MeV. 

The intermediate situation {cr[^'> ~ as) can be interpreted as a statistical 
equilibrium between capture and stripping. For a given combination {Zion, 
Zmedia, Eion) ouc cau define a probability p{q) for the ion to have a charge q 
(or, in other words, to travel with Z — q electrons). 



2.2 Stripping materials 



Considering equation 2 and the fact that we have no possibility to directly 
determine the charge state of an ion, it becomes clear that the mass identi- 
fication of fragments requires, as far as possible, those fragments to be fully 
stripped in the spectrometer. To optimize the proportion of bare ions, it is 
convenient to insert layers of well-chosen materials in the device. The choice 
of a stripper foil results in a compromise between its stripping efficiency and 
its impact on the beam quality, in terms of angular and energy straggling, and 
in terms of induced secondary reactions in the stripper. 




Projectile energy (AMeV) Projectile energy (AMeV) 




Projectile energy (AMeV) Projectile energy (AMeV) 



Fig. 2. Bare ion proportion (left scale, continuous lines) and nuclear reaction rate 
(right scale, dotted lines) in various stripper foils (see legends) for several heavy 
beams: Mo, Xe, Pb and U, as a function of the incident energy. The thickness chosen 
for the stripper corresponds to its charge equilibrium thickness at each energy. 
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Figure 2 displays the stripping efficiency, evaluated from the fraction of bare 
ions, as function of the projectile energy. Calculations were made using the 
code GLOBAL [5]. The percentage of nuclear reaction in stripper is also dis- 
played (reaction cross sections were estimated using the Karol parametriza- 
tion [6]). In those calculations, to simplify the discussion, the stripper thick- 
ness has always been chosen as the so-called charge- equilibrium thickness. 
This quantity corresponds to the thickness of matter beyond which the pro- 
portion of the ionic charge states is fully independent of the initial charge-state 
distribution. 

At high energy, any high- or medium-Z element will provide excellent results. 
Nevertheless, even at lA GeV (the highest energy available at the GSI syn- 
chrotron SIS), results are not perfect for very heavy elements: fully stripped 
ions represent only 95% of Pb ions, and no more than 88% for U. 

With decreasing energy, the efficiency of strippers is subject to large variations. 
Nb is the most efficient stripper for beam energies higher than 600A MeV, 
while Ti gives slightly better results below this value. Nevertheless, the equi- 
librium charge state is reached in a thicker foil in Ti than in Nb, which results 
in a higher nuclear-reaction rate. Let us point out that the energy regime in 
which the situation can be considered as favorable varies strongly according to 
the considered ion beam. The 90% proportion of fully stripped ions, which can 
be considered as a reasonable limit, is hardly reached for a lA GeV U beam, 
while it allows experiments with a Xe beam at energies as low as 2 00 A MeV. 



2.3 Impact of charge states on the separator selectivity 

The thick degrader technique has been developed in order to allow the produc- 
tion of a high-purity secondary beam. This technique consists in the inclusion 
of a thick layer of matter between two magnetic selection devices. The first se- 
lection basically operates on the A/q ratio of the ions. Those ions loose a large 
part of their energy in the degrader, which basically depends on the square 
of their nuclear charge. The second magnetic selection can then be roughly 
understood as a Z selection. 

The selectivity of such a separation device is presented in figure 3. Galcula- 
tions have been performed in order to simulate the transmission in the FRS 
for fragments produced in the Pb-l-p reaction at 500A MeV. We only took 
into account the fragmentation-evaporation products and did not consider 
the production cross sections of the reaction products (these cross sections are 
expected to drop very rapidly for mass losses beyond 40 and completely vanish 
beyond 60). From the hundreds of different fragments that we considered, less 
than 20 are transmitted up to the end of the separator. The use of shts in the 
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Fig. 3. Simulation of ions transmitted by the FRS spectrometer. The beam is Pb 

with an energy of 500 A MeV. A thick degrader of thickness 2000 mg.cm~'^ Al is set 
between the two magnetic sections. Fragments transmitted in the first part of the 
spectrometer are plotted in grey, while those also transmitted in the second part are 
plotted in black. Only fully stripped ions are considered on the left figure, while all 
combinations with or 1 electron in the first part and 0, 1, 2 or 3 electrons in the 
second part are considered on the right figure. 



middle of each magnetic sections would even reduce this number. 

Things are less favorable when ionic charge states are taken into account. 
The various combinations of charge states allow for a much larger population 
of fragments to be transmitted, as can be seen on the right part of figure 3. 
The purity of a secondary beam is strongly reduced. Consequently, production 
measurements require more time to achieve satisfactory statistics as a large 
part of the transmitted fragments may not be of interest. 

It is important to notice that, due to the bell shape of the isotopic distribu- 
tions of evaporation or fission residues, the situation is drastically different if 
one considers proton-rich or neutron-rich isotopes. If the magnetic device is 
set to select fully stripped fragments with a given A/Z ratio, H-like fragments 
with atomic number Z and mass A — A/Z will also be transmitted, as well as 
He-like fragments of mass A — lAjZ ^ and so on. The proportion of unwanted 
ions therefore directly depends, not only on the charge-state probability, but 
also on the production cross sections for the lighter fragments. Let us consider 
fragments of a given element. In the case of neutron-rich isotopes, cross sec- 
tions for lighter isotopes are considerably larger, and therefore the proportion 
of contaminants is expected to be very important. On the other hand, when 
considering proton-rich fragments, the cross sections for the production of the 
lightest isotopes are very small, and therefore the problem of charge states 
may be somehow disregarded. 
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3 Z identification 



3.1 Set-up design 



750A MeV, 1 MUSIC (1 bar) 



750A MeV. 2 MUSIC (1 bar) 




■ Z=82 

■ Z=81 

■ Z=80 
Sum 
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350A MeV, 4 MUSIC (2 bar) 




Fig. 4. Simulation of the energy loss of tfiree fragments (Hg, Tl, Ph) in MUSIC 
chambers [7j (60 cm long each) filled with 1 or 2 bar PIO gas. The fragments 
were assumed to enter the detectors fully stripped. Incoming energy of fragments 
is 750 A MeV (upper part) or 350A MeV (lower part). 

Calculations presented in figure 4 summarize the different possible behaviors 
for stripped ions entering an ionization chamber. The setup in our simulation 
was intended to simulate the MUSIC chambers [7] used at the FRS. In normal 
conditions, these chambers are filled with PIO, a gas mixture with 90% of Ar 
and 10% of CH4, at atmospheric pressure (this corresponds to a gas thickness 
of 100 mg.cm"^). As only 2/3 of these chamber length is instrumented, the 
ions were propagated along all the gas thickness, but only 2/3 of the thickness 
was used to calculate the energy collected after deposition in the gas. 

The electron stripping and/or capture probability was computed from the 
cross sections as calculated by GLOBAL [5] for each small step of the fragment 
propagation through the gas, and the energy loss calculated according to the 
new ionic charge state of the fragment. Effects such as the escape probability 
of the S electrons or the drift of the electrons through the chambers were not 
taken into account. The contribution of effects other than the variations of the 
ionic charge state were roughly taken into account as a fixed 0.3% contribution 
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to the total straggling. 



On the upper part of the figure, the incoming ions have an energy of 750A MeV 
and pass through a limited gas thickness, which corresponds to a single ion- 
ization chamber filled with 1 bar gas. Their mean free path related to electron 
capture is much larger than both the gas thickness and the mean free path 
related to electron stripping. The ions remain stripped along most of their 
way through the gas, and the achieved resolution in energy loss is excellent. 
The small tails associated with each peak correspond to the ions that kept 
an electron for some time in the gas. Of course, adding a second chamber 
increases the resolving power of the system. 

If the incident energy of the fragments is lowered to 350A MeV, the probabil- 
ity of charge exchanges during the crossing of the fragment through the gas 
dramatically increases. Therefore, the energy-loss distribution becomes wider 
as a result of the fluctuations of . leading to an unacceptable decrease of the 
resolution of the ionization chamber. If one increases the total gas thickness 
(up to 700 mg.cm"^ on the figure) the relative statistical fluctuation in is 
reduced, leading to a better definition of the total energy loss and therefore 
a better resolution on . According to this simulation, a resolution better 
than 10% may be expected. Furthermore, complete independency regarding 
the incoming charge state is achieved. 



3.2 Confrontation with experiment 

According to simulations, the use of ionization chambers to identify heavy, low- 
energy fragments seems possible if one strongly increases the amount of gas (or 
the number of chambers) with respect to high-energy setups. This technique 
has been used successfully during the Pb-l-p at 500A MeV experiment at GSI. 
The setup consisted of 4 MUSIC chambers [7] filled with 2 bar PIO, for a 
total gas thickness of the order of 800 mg.cm"^, from which roughly 2/3 were 
instrumented. Due to the presence of a thick degrader at the intermediate 
focal plane, the energy of fragments in the second part of the spectrometer 
was of the order of 300A MeV. The resolution obtained in nuclear charge was 
better than 0.9% (or 0.7/Z) for very heavy fragments and better than 0.6% 
(or 0.4/Z) for fragments around nuclear charge 70. The resolutions given here 
are FWHM. We will keep this convention throughout the paper, unless the 
contrary is specified. 

Examples of resolutions obtained are presented on figure 5 for fragments close 
to (left part of the figure) and far from (right part) the lead projectile. This 
resolution obtained by the combination of the signals of the 4 chambers can 
be compared with the resolution obtained with a single chamber, which can 
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z z 

Fig. 5. MUSIC spectra obtained in the Pb+p experiment at 500 A MeV for FRS 

settings centered on ^^"^Tl (left) and Yb (right). 

be deduced from figure 7. To obtain these results, a correction taking into 
account the velocity of fragments had to be applied. 

But the large gas thickness had three drawbacks. 

(1) The energy loss is large enough so that a dependency regarding to the 
mass of the fragment was observed. This is a trivial expectation issued 
from the Bethe formula, but the consequence is that, for a given Z, a 
systematic increase of the MUSIC signal was observed with decreasing 
masses. The experimental evidence for this effect is presented in figure 6. 
One can clearly see that this effect is strong enough to lead to a misiden- 
tification of the fragments by one charge unit, if not taken into account. 

(2) The large amount of matter increases the parasitic nuclear-reaction rate 
in the gas. Products from these reactions can be rejected during the off- 
line analysis by comparing the signals observed in the different MUSIC 
chambers. An example of such a correlation is presented in figure 7. 

(3) The large quantity of energy deposited in the gas seemed to strongly 
alter the shape of the relation between the signal of the MUSIC and the 
position of the fragments in the chambers, as well as the resolution (see 
section 3.3). 

3.3 Effect of large energy deposition in ionization chambers 

In ordinary operating conditions, a reduction of the signal is observed with 

increasing distance between the point of passage of the fragment trajectory 
and the anodes increases. This dependency is due to the recombination of 
some of the drift electrons in the gas along their path toward the anodes. 
Those recombinations may occur with PIO molecules, but are more likely to 
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Fig. 6. Mean MUSIC signal for several isotopes in the Pb+p experiment at 
500 A MeV. Each abscissa corresponds to a given setting of the magnets of the 
FRS. In first approximation this defines a cut on the mass-over-charge ratio of the 
fragments. Considering the limited acceptance in magnetic rigidity of the FRS, each 
point corresponds to an average over 3 to 5 neighboring masses. 




1400 1600 1800 

Signal of MUSIC 1 [a.u.] 

Fig. 7. Correlation of signals from the first and the second MUSIC chambers in a 
setting centered on ^^^Yb in the Pb+p experiment at 500A MeV. 

occur with impurities such as O2 or H2O. The left part of figure 8 presents an 
example of this effect. 

If the counting rate is increased up to some 2 kHz, a drastic change of the 
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Fig. 8. Sum of the energy deposition in the 4 MUSIC chambers as a function of 
the distance from the anodes (arbitrary units with zero set in the middle of the 
chambers) in the Pb+p experiment. Count rate is of the order of 300 Hz (left) and 
2 kHz (right). 



Time [u. a.] 



Time [u. a.] 



Time [u. a.] 



Time [u. a.] 



Fig. 9. Sum of the signals of the 4 MUSIC chambers as a function of the time 
duration of the SIS spills, integrated along different parts of the axis transversal to 
the trajectory of the fragments. 8 regions of equal length have been defined on this 
axis. From left to right the figures correspond to the first, third, fifth and seventh 
region, the first being the closest to the anodes. The counting rate is in the order of 
2kHz. 

relation between the signal and the drift distance is observed (right part of 
the figure 8). Also, compared to the figure on the left, if the resolution is only 
slightly altered in the region close to the anodes, no identification is anymore 
possible for trajectories far from the anode. 



Further illustration of the relation between count rate and resolution in the 
case of high event counting rates is presented on figure 9, in which the MUSIC 



13 



signal is plotted as a function of the time duration of a beam spill of the SIS 
synchrotron. For a section of the MUSIC close to the anodes, the signal slightly 
reduces with increasing beam intensity and then gets back to its nominal value 
at the end of the spill. This limited effect probably is the consequence of the 
appearance of some screening effect. As one gets further from the anode, the 
relation inverts: increasing intensity leads to increasing signal, up to a point 
that leads to a misidentification by 2 nuclear charges. As this effect depends 
on many parameters it could not be corrected during the analysis, resulting 
in a loss of some data when the identification of the ions becomes impossible. 

As this variation of signal depends on position, electronics can not be responsi- 
ble for it. This effect is hkely to come from the saturation of electron catching 
sites of the gas (mainly impurities) for high electron densities corresponding 
to a large number of highly ionizing nuclei. 

Such short time variations on MUSIC signals have already been observed in 
the past without satisfactory explanation. They were tentatively ascribed to 
pressure or temperature drifts, but no clear correlation with such quantities 
could be established. We have shown here that this behavior is hkely to be 
related to the electron-attachment saturation under high counting rate, or in 
other words to the amount of ionization in the chambers. 



3.4 An original method for heavy fragment identification in Z at low energy 



We have described a new way to operate the ionization chambers. We have 
shown that such a set-up allows to identify in Z heavy ions as heavy as lead 
at energies as low as 300A MeV. We are confident that this method would be 
efficient for lower energies. 

The limitations of this set-up have also been stressed. First, the large amount 
of gas induces variations of the energy loss with the mass of the fragment. This 
effect does not hamper the efficiency of the method but must be taken into 
account during the analysis process. Second, as the density of ionization in 
the gas largely increases with decreasing energy, the response of the ionization 
chamber can be modified according to the counting rate. This effect can be 
strong enough to prevent any identification of the fragments. Therefore, it may 
be necessary to limit the counting rate. 
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4 A thick degrader as a passive measurement device 



As written in equation 2, the mass identification of a fragment requires the 
knowledge of its ionic charge state. To the best of our knowledge, no apparatus 

can provide a direct, unambiguous measurement of the ionic charge state of 
an ion in the energy range considered in this paper. Nevertheless, it may be 
deduced from a fine evaluation of the energy loss in a thick layer of matter, 
namely the degrader set between the magnetic sections of the FRS. 

In this section, we will consider that the fragments are correctly identified 
in Z, regardless of their charge state, according for example to the use of 
the method presented in the previous section. For the sake of simplicity, we 
will often call "degrader" the whole matter inserted in the beam line at 5*2, 
namely, the degrader plus the plastic scintillator. They were of course treated 
separately in the simulation. 



4-1 Energy loss in the degrader 



The energy loss in matter of an ion of mass A and kinetic energy E is given 
by: 

AE = AA7 (3) 



For calculations related to magnetic rigidities, it is more convenient to write 
equation 3 as a function of A{Pj) instead of A7. This leads to: 

AE ~ AA((3-f) <(3> (4) 



where < /3 > is the mean value of (5 in the considered layer of matter. In the 

energy domain wc are discussing in this paper, the variations of this quantity 
are small and one can therefore treat it as a constant for the purpose of the 
demonstration. 

The slowing down of the fragments and the change of charge state (if any) 
induce a variation of the magnetic rigidity between magnetic sections. Using 
the discussion that followed equation 2, this can be written as: 



A(Sp) = -(/37)i--(/37)2 (5) 

= - ((/37)i - {h)2) - f-(/?7)i - -(/97)i) (6) 
0.2 \Q2 Qi J 
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Schematically the first A/q2 term (which we will note as A(i?p)(^)) corre- 
sponds to the loss of energy, while the second term (which we will note as 
A(i?p)('^)) is related to the charge-state changing. 

According to equation 4 one can write: 

A(5p)(-) = ^1 (7) 



According to the discussion in section 2, AE can be written: 

AE = qlfff{P)g{Xs,) (8) 

where ggjj is the mean value of the square of the atomic charge of the fragment 
in the degrader (as defined in chapter 2), f{f3) is the velocity dependence from 
the Bethe formula, and g{Xs2) describes the position dependence due to the 
variations of the degrader thickness. 

According to its definition, g^^^ is completely independent of the ionic charge 
state in both parts of the spectrometer. Furthermore, it is both unambiguously 
related to and very close to it. We can then replace it by for this 
calculation. We can now write A{Bp)^^'> as: 

= ^.Z + e. (9) 

?2 F q2 



Here 62 stands for the number of electrons in the second section of the spec- 
trometer. The second part of the development of A[Bp) writes: 

A(Sp)W = {Bp),.^^^ (10) 

Q.2 



To summarize, we have established that, besides the initial velocity of the 
fragment, the variation of magnetic rigidity depends on three physical com- 
ponents: 

• the nuclear charge of the fragment; 

• the charge state of the fragment in the second part of the spectrometer; 

• the change of charge state of the fragment between the two sections of the 
spectrometer. 

In the next section (4.2) we will present the possible applications of the mea- 
surement of magnetic- rigidity variation and the associated requirements on the 
energy-loss resolution. In section 4.3 we will discuss the limitations from the 
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characteristics of the devices and from physics. In section 4.4 we will present 
and discuss the results obtained in the 500A MeV experiment. 



4-2 Applications of the measurement of the energy loss in the degrader 

Among the three effects listed in the previous section, the charge-state chang- 
ing is always dominant as soon as a capture or a ionization takes place in the 
degrader, even with a very thick degrader for which the energy loss is large. 
This allows to detect and identify unambiguously any change of charge state 
as will be shown later in figure 12. The related magnetic-rigidity variation 
does not depend on the thickness of matter set, and therefore information 
on a change of charge state is accessible whatever the setup or the energy of 
the incoming fragment. During the data analysis, the fragments can be sorted 
according to their change of charge state, a quantity that we will write as Ag. 

Besides this property, and as expected from equation 9, a large thickness of 
matter induces a linear dependence of the variation of magnetic rigidity with 
Z and q. Separating peaks associated to neighboring charges with a reasonable 
efficiency requires the FWHM of each peak to be (roughly) smaller than 0.7/Z. 
For example in the case of lead this corresponds to a resolution of 0.9%. 

In high-energy experiments, the selection of the charge states of the fragment 
is nearly entirely obtained using the Ag value. Indeed, although the propor- 
tion of fully stripped ions (that we can write (0,0) according to the number of 
electrons they carry in each part of the spectrometer) is only of the order of 
80%, the Ag estimation allows the rejection of most of the charge-state com- 
binations like (0,1), (0,2) or (1,0). The only charge-state combinations which 
can not be disentangled from (0,0) are (1,1), (2,2), and so on. But in the case 
of the fragmentation of a lA GeV Pb beam for example, these charge states 
represent only 2% of the Pb fragments that have a Ag equal to zero. This 
very low proportion is due to the combination of the two strippers (one at the 
entrance of each magnetic section) which is weighted by the product of the 
two probabilities. These fragments might be disregarded during the identifi- 
cation process and receive an appropriate correction during the estimation of 
production rates. In this case it is therefore convenient to use the degrader as 
an additional measurement of the Z value, which can be combined to the one 
obtained from the MUSICs and increase the overall Z resolution. This method 
has already been successfully used in previous spallation experiments at the 
FRS [8,9,10,11,12]. 

In the case of a low-energy experiment, the fraction of Ag = fragments can 
be as low as 40% for fragments produced from a 500A MeV Pb beam. Among 
those Ag = fragments, roughly 40% are hydrogen-hke ions. Therefore, they 
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should be explicitly identified and taken into consideration. If one relies on the 
MUSICs for the measurement of the nuclear charge, according to the previous 
calculations, a precise estimation of the energy loss in the degrader may be 
used to estimate the ionic charge state of the fragments. 



4-3 Resolution limitations imposed by physics and by detectors 

Let us first stress what exactly means the word resolution. The value we 
are considering is the variation of the magnetic rigidity. It is important to 
distinguish the two sources of the so-called uncertainty on this quantity. The 
first and most obvious one is the uncertainty on the measurement. The second 
one, by far the largest as we will see, is not an uncertainty in the usual sense of 
the word: it is the fluctuation in energy induced by the path through a thick 
layer of matter. Whatever the precision of the measurement of the magnetic 
rigidity or time of flight, there is no way to correct the ABp width induced by 
those fluctuations. Therefore, to identify charge states, it is the sum of those 
two terms which must fulflU the 0.7 /Z condition presented above. 

For convenience, one can express the straggling in energy as a fluctuation of 
the variation of the magnetic rigidity: 



This relation is discussed in details in appendix A. In the following we will 
express all fluctuations and uncertainties as relative uncertainties with respect 
to the variation of magnetic rigidity. 

We will now discuss separately the contributions to the total ABp width of 
the energy straggling and of the magnetic- rigidity measurement. 



4.3.1 Measurement of the magnetic rigidity and associated uncertainty 

The resolution that can be achieved on the measurement of the change of 
magnetic rigidity obviously depends on the resolution in magnetic rigidity, 
which itself depends on the precision of the magnets and from the position- 
measurement devices, namely scintillators. 

In a given magnetic setting between the focal planes A and B, characterized 
by a magnetic rigidity (-Bp)o for centered ions, the magnetic rigidity of an ion 
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is deduced from its position at the focal planes (respectively xa and xb) as: 
(Bp) = {Bp)o. (l + ^ - M^) (12) 

with D being the dispersion of the magnetic section and M its magnification. 

An optimum result is obtained only if the magnetic optics exhibits no aber- 
ration and if the positions of the fragments are measured at the focal planes 
of each magnetic section. We will assume that the first condition is fulfilled, 
although there are indications that some limited magnetic aberrations may 
exist close to the final focal plane. 

For the second assumption, the situation differs according to the focal plane 
one considers. In the first part of the FRS, the initial point (the target) is 
unambiguously defined, and the scintillator is set at the intermediate focal 
plane, at a fixed position. The magnification term in equation 12 vanishes and 
the uncertainty is simply obtained from the uncertainty on magnetic fields 
(known to be less than 10~^) and on the scintillator response (estimated to be 
at most 3 mm): 

(Sp)i Bi D ^ ' 



The position of the final focal plane varies according to the chosen magnetic 
optics, and the position of the corresponding scintillator depends on the avail- 
able room left by all other detectors located at this region. Therefore, the 
two points may not coincide. If the distance between the scintillator and the 
focal plane is d and the angular distribution of a fragment has a width ae, a 
dispersion is introduced on the measurement of the X4 position: 

6x4 — ag.d (14) 



The contribution of this term to the uncertainty on the measured magnetic 
rigidity adds to the intrinsic uncertainty of the measurement: 

If the trajectory angle of the fragments is measured, the position of fragments 
at the final focal plane can be reconstructed. In this case, the enlargement on 
the position resolution due to the incoming angle can therefore be corrected. 
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4-3.2 Energy straggling 



Two effects are responsibfe for tlie fluctuations of tlie energy loss of a fragment 
passing through a solid [14,15]. First, the energy loss is due to successive 
collisions, each removing a variable amount of energy from the incoming ion. 
The variance can be estimated with a reasonable precision by: 



where is the nuclear charge of the material, Zf the nuclear charge of the 
fragment passing through, and Na the number of atoms of the material. This 
term will decrease the resolution as the thickness of the material increases 
and as the energy of the fragment increases. In the Bohr approximation x 
reduces to 1-/3^/2. The Lindhard-Sorensen (LS) theory [13] has extended the 
validity of this equation to relativistic energies. In this theory x is a much more 
complicated expression but the agreement with experiment in the relativistic 
regime is strongly improved. For our subsequent calculations of energy losses 
and straggling, we used this theory as implemented in the code ATIMA [16]. 

The second phenomenon inducing straggling is the fluctuation of the ionic 
charge state, already discussed in section 3. These fluctuations are of statistical 
nature, therefore the energy loss resolution increases with the thickness of 
the solid. On the other hand, fluctuations increase with decreasing energy as 
more and more charge states can be populated. The calculations related to 
this phenomenon were performed using the ionic capture and stripping cross 
sections from the code GLOBAL [5]. 

We have performed calculations of these two effects for ions of ^o^Pb (left 
part of figure 10) and ^'''^Hf (right part). For each energy we calculated the 
energy loss of a ^osp]-, beam in a hquid hydrogen target used in all spallation 
experiments performed at the FRS (see chapter 1), but did not consider any 
straggling or angular broadening at this point. Any energy straggling in the 
target is supposed to be washed out by the achromatism of the system. This 
is not the case of the angular distribution induced by the reaction, but we 
disregarded it in order to focus on the effects of the matter set at S2. The 
thickness of the degrader was then calculated for each energy as it is done 
during experiments: the total thickness of matter at the intermediate focal 
plane must correspond to half of the range of the beam in aluminium (this 
value represents an optimum between the energy loss and the probability for 
parasitic nuclear reactions in these layers). We explicitly took into account a 
plastic degrader made of 3 mm of BC-400. 

Results of these calculations are expressed in terms of FWHM for the relative 
magnetic rigidity variation. At lA GeV the charge-state ffuctuations are com- 




(16) 
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Fig. 10. Calculations of the straggling on the energy loss (FWHM) due to collision 
straggling (circles) and charge state fluctuations (squares) for "^^^Ph (left plot) and 
^"^^ Hf (right plot) , as the function of the energy of the beam. The continuous line rep- 
resents the quadratic sum of these two components. Effect of the angular straggling 
is presented in the hypothesis of a distance of 1 m between the position measurement 
at the exit of the spectrometer and the position of the magnetic focal plane. See text 
for details of the simulation. 

parable to the effect of tlie collision straggling in the case of Pb, while, due 
to the quadratic summation, they are nearly negligible in the case of Hf. The 
total FWHM for Pb ions is 0.49%, and 0.44% for Hf. 

With decreasing energy, the collision straggling keeps a nearly constant value, 
while the change state changes intervene more and more. At 500A MeV, the 
resolution on energy loss for lead is not better than 0.84%, the collision strag- 
gling accounting for only 20% of this value. 

Unexpectedly, in the case of Hf, the worsening of the resolution with decreasing 
energy is stronger than in Pb and does not reach a maximum as it does for 
Pb. This is probably related to the number of collisions in the materials, 
which increases faster for Pb than for Hf with decreasing energy, therefore 
reducing the statistical fluctuations on the charge state. Furthermore, around 
250A MeV, the cross section for capture to the s atomic shell of Pb ions is 
already very large while the capture cross section to the /-shell remains small: 
this strongly favors the q = Z — 2 charge state and thus furthers reduces the 
fluctuations. In the case of Hf this situation appears only at an energy lower 
than 200 MeV. 

We did not extend the calculations to lower energies as the energy loss in the 
target would have become too large to render an experiment feasible. Reducing 
the hydrogen thickness would not counteract this problem in a relevant way 
as it would also lead to an increase of the relative contribution of the target 
windows (which cannot be reduced) to the production rates. 

We also calculated the effect of the angular straggling induced by the materials 
set at the intermediate focal plane. In chapter 4.3.1 we have already discussed 
this effect: in an achromatic system like the FRS, it varies linearly with the 
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Fig. 11. Beam width in percent of deviation of magnetic rigidity, as deduced from 
the positions measured at S2 and S4 in the Pb+p experiments at 500A MeV (left) 
and lA GeV (right). Squares correspond to experimental points, and a Gaussian fit 
has been performed on each spectrum (solid curve). 

distance between the real magnetic focal plane and the point of the beam line 
at which the position measurement is performed (it is therefore canceled if the 
position is measured at the correct point of the beam line). The values plotted 
in figure 10 correspond to a distance of 1 m between the expected and the 
actual measurement points. As one can see on the figure such a misplacement 
has nearly no consequence on the obtained energy loss resolution. Of course a 
larger distance could play a significant role. 

4.3.3 Experimental beam width at the final focal plane 

As we have already discussed, in an achromatic system, the position of an ion 
at the final focal plane does not depend of its initial velocity, nor of its initial 
angle. Going through the degrader induces velocity fluctuations that can not 
be compensated by the magnetic optics. Therefore, the beam is not expected to 
be perfectly focused at 5*4. The previous calculations on the energy straggling 
can be converted in a minimal width expected for the beam at 5*4: the results 
are 0.27% at 500A MeV and 0.16% at lA GeV. The experimental values are 
presented in figure 11. Beam widths at S2 are also plotted for comparison ; 
let us recall that the dispersion observed at the intermediate focal plane is 
expected to be fully corrected by the achromaticity of the system. 

In the case of the lA GeV experiment, the width observed at S2 corresponds 
to a magnetic-rigidity straggling of 0.06% (FWHM relative to Bpi). An en- 
largement is observed at 6*4, leading to a magnetic-rigidity straggling of 0.2% 
(FWHM relative to Bp2, thus 0.54% relative to ABp). This is in good agree- 
ment with the result of the calculation presented above (0.49% FWHM relative 
to ABp). 

In the case of the 500A MeV experiment, the width observed at S2 corre- 
sponds to a magnetic-rigidity straggling of 0.16% (FWHM relative to Bpi). A 
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very large value is observed at 5*4. corresponding to a magnetic rigidity strag- 
gling of 0.4% (FWHM relative to Bps, thus 1.24% relative to ABp). This is 
much higher than the result of the calculation presented above (0.26% FWHM 
relative to Bp2)- Let us consider the different possible reasons for this large 
position spread. 

• The absence of correction for the angle of the ions. In this discussion we 
only consider the beam and therefore this possibility can be ruled out. Even 
if the scintillator is far from the magnetic focal plane (2 m can be considered 
as an upper boundary) , the position broadening induced by the angle is not 
expected to be a major contributor to the overall broadening (the angular 
aperture of the beam has been estimated by our calculations to be of the 
order of 2 mrad at this point of the beam line). 

• The energy straggling and the charge state fluctuations. This obviously 
plays a role: we have shown in section 4.3.2, that precise calculations lead 
to a value of 0.26% (FWHM relative to Bp2) for this effect. This induces an 
important broadening of the position spectra (which is not compensated by 
the achromaticity of the system) but is not sufficient to explain the observed 
value. 

• A failure in the design of the degrader leading to a chromaticity of the 
system. If one rules out an error in the calculation of the degrader proflle, 
the hypothesis of surface defects or inhomogeneities remains. 

Prom the difference between the calculated widths at S4 at the two considered 
energies, we can postulate a typical size of the inhomogeneities. We could re- 
produce the observed widths at both energies by considering that the degrader 
surface is inhomogeneous with very small defects: 27 fim. This value is to be 
compared to the degrader thicknesses: roughly 6.5 mm in the 500A MeV ex- 
periment and 15 mm at lA GeV. The relative straggling in magnetic rigidity 
difference induced solely by those inhomogeneities would be 1% at 500A MeV, 
and only 0.3% at lA GeV. 

We consider the presence of such inhomogeneities plausible as their effect 
would be small enough not to affect experiments conducted at higher ener- 
gies, which necessarily include thicker degraders. Therefore they were never 
considered as potentially harmful for experiments. 



4-4 Obtained resolution on the variation of magnetic rigidity in high- and 
low-energy experiments 

The use of a thick degrader as a passive measurement device has already 
proved to be successful in previous experiments conducted at the FRS at 
energies of or close to lA GeV. We have already emphasized that its use 
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Fig. 12. Correlation between the variation of magnetic rigidity (normalized to the 
nuclear charge for fully stripped fragments) and the Z value obtained from MUSICs 
in the Pb+p at 500 A MeV experiment. The central band corresponds to fragments 
that did not gain or loose any electron while passing through the degrader. The upper 
band corresponds to fragments that lost an electron, the lower band to fragments that 
gained one electron. 

in the spallation of lead at 500A MeV was of particular importance as the 
number of non-stripped ions was expected to be very high (up to 60% after 
the degrader). 

At higher energies, the influence of degrader inhomogeneities become negli- 
gible, and the charge-state straggling is strongly reduced. Our calculations 
indicate that a 0.5% resolution should be obtained, which allows for a satis- 
factory Z or q identification. This conclusion is in good agreement with results 
from previous experiments. 

From figure 12, it is obvious that during the 500A MeV experiment the reso- 
lution required to separate charge states was not achieved: no separation can 
be seen between fully stripped and H-like fragments. They appear as a bump, 
clearly visible for each spot in the upper part of the most intense line, but 
impossible to disentangle from the fully stripped fragments. This could be 
expected from the discussion conducted in the previous section: the position 
straggling observed at 5*4 leads to an uncertainty of 0.4% on the magnetic 
rigidity. As the difference of magnetic rigidity is only 31% of the magnetic 
rigidity in the second part of the spectrometer, the uncertainty on ABp is 
expected to be close to 1.25%, which is very close to the observed resolution. 
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Such a value prevent identification from the various charge states of the frag- 
ments. In the discussion above we only considered the beam, which allowed us 
to neglect the influence of the angle on the precision of the magnetic-rigidity 
variation. This is not true anymore when one looks at the fragments, and a 
non-optimal position of the scintillator at 84^ may have further worsen the 
experimental resolution. We could not put in evidence any effect of the angle 
in the case of the 500A MeV experiment as a failure of a device prevented 
the measurement of the angle of the fragments at 5*4. This does not affect our 
global conclusion on the reasons why the reduced resolution was obtained for 
the energy loss measurement in the degrader. 

Based on our simulations and the experimental results, we postulate that 
defects on the degrader surface are the main reason for this very poor resolu- 
tion, and that an improved degrader (with a surface roughness smaller than 
10 /im) may allow this method to be successfully used at an energy as low as 
500A MeV for very heavy beams such as lead. 



5 Evaluation of fragment masses 



The determination of the mass of each fragment requires two steps, which will 
be described in this section. First, the A/q ratio is deduced from the magnetic 
rigidity and time of flight in the second section. Then, if the use of the degrader 
has been successful in determining the ionic charge state of the fragments, or 
if the proportion of non-stripped ions is low enough to neglect them, getting 
the mass of each fragment is straightforward. 

Alternatively, if the charge state is not known, it is convenient to determine a 
"most-hkely" mass for each fragment. 



5.1 A/q ratio 



The A/q ratio is deduced from the magnetic rigidity and from the time of flight 
in the second section of the spectrometer. Various corrections of dependencies 
from the positions have to be applied: effect of the position of the fragment 
at S2 and 5*4 on the length of the trajectory, non-linear effects in the time 
response of the scintillators. This leads to a parametrization defined as: 

A/q = fix,, X4, ToF) = \^^'''':l^f' , (17) 
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Li and T, arc free parameters. They are assumed to stay constant as long 
as magnetic optics is not changed. Their vahies can be adjusted considering 
a given setting for which masses are affected by counting, starting from the 
centered fragment. It is usually convenient to use a magnetic setting corre- 
sponding to fragments close to the beam, as identification is in this case quite 
straightforward because of the reduced velocity spread of those fragments. In 
contrast, in low-energy experiments the high proportion of charge states sug- 
gests a setting centered on very light isotopes. The calculated parameters can 
be checked by re-calculating the mass number of the ion beam. In the case of 
the 500A MeV experiment, the mass of ^°^P6 was obtained with a precision 
of 5.10-1 



5.2 A/q resolution and effect of the polluting charge states 

The resolution required to separate masses with a peak-over-valley ratio of 
10 is 0.1%. The FRS layout and detectors have been designed to fulfill this 
requirement. In this respect, the most difficult situation is expected to be 
the measurement of heavy fragments at high energy: in this case the time of 
flight is short and the difference in magnetic rigidity between two neighboring 
masses is minimum. In this respect, a low-energy experiment is expected to 
be a better the first constraint is softened. 

Nevertheless, the presence of polluting charge states can strongly hamper the 
resolution obtained on the A/q ratio. For fragments produced from a heavy, 
stable nucleus such as lead, the A/q ratio is close to 2.5. Therefore the A/q 
peaks for fragments with an odd number of electrons will fall in-between peaks 
associated to fragments with an even number of electrons, leading to an ap- 
parent loss of resolution. This has of course no consequence if the ionic charge 
state of each fragment is known, as the exact mass can be obtained. On the 
other hand, if the ionic charge state is not known, one can assign its most- 
likely mass to each fragment. Also, according to the previous discussion, the 
fragments of mass A are mixed with fragments A — 2 and A — 3. We will show 
that one can use the correlation between the A/q ratio and the ABp values 
to reduce the number of actual masses that are mixed together. 

An example of such a correlation is presented in figure 13, for heavy and light 
Hg isotopes produced in the 500A MeV Pb-l-p experiment (only fragments with 
Aq = are represented). This figure also underlines the lack of resolution in 
ABp obtained in this experiment. The large bands seen for heavy isotopes are 
a mix of the fully stripped ions (low part of bands) and H-like ions (upper 
part). In the case of light isotopes, very few H-like ions are produced and 
transmitted, and therefore only the lower part of the bands remains. 
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Fig. 13. Correlation between the differences of magnetic rigidity and the calculated 
A/q ratio for Hg ions. Left figure corresponds to a FRS setting centered on ^^^Hg, 
right part to a setting centered on ^^^Hg. 

As no direct separation of charge states is possible, it is advantageous to 
get the A/q value, not from a direct projection on the A/q axis, but from 
a projection along the direction of the correlation bands. Not only does it 
improve the apparent resolution, creating clearly separated peaks (in the case 
of the 500A MeV experiment, the resolution becomes close to 0.1% for all 
settings, while it could reach 0.4% in the most unfavorable cases). It also 
ensures that only masses A — 3 contribute to the pollution of isotopes of mass 
A. The unavoidable weakness of this method is that the group of fragments 
identified as being of mass A and atomic charge go contains a variable fraction 
of {A — 3,q = qo — 1) fragments, which will have to be subtracted in order to 
deduce the production rates. 

At this step unambiguous identification in mass on an event-by-event basis is 
of course not possible anymore. Nevertheless it is convenient to assign a mass 
to each fragment. To calculate masses, we assume that fragments gaining zero 
or more electrons when passing through the degrader were fully stripped in 
the first part of the FRS, while fragments loosing one or more electrons are 
fully stripped in the second part of the FRS. In other terms, for a given value 
of Ag, each fragment is assumed to be in the most probable charge state in the 
two parts of the FRS. In the worst case (heavy fragments for which Aq = 0) 
this hypothesis is valid for only 60% of fragments. 



6 Production rates for ambiguously identified fragments 



Even if fragments are not unambiguously identified in mass on an event-by- 
event basis, we will show that it remains possible to extract isotopic production 
rates. This requires to evaluate and correct for the fraction of fragments that 
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was not correctly identified because of tlie unresolved ionic charge states. In 
this chapter we will present a method to achieve those corrections. 



6.1 Kinematic spectra 



First of all it is necessary to take into account the limited acceptance of the 
FRS. Its acceptance is limited to ±1.5% in magnetic rigidity (or in longitudi- 
nal momentum) for each given ion. This means that, if the broadening of the 
fragment velocity goes beyond this value, only part of the momentum distri- 
bution will be transmitted in a given setting of magnets. Calculation of the 
production cross sections therefore requires the reconstruction of the complete 
kinematic distribution of each fragment, which may require data from different 
magnetic settings. 

The momentum from each fragment in the first part of the spectrometer can 
be calculated by applying the relation (2) to the first part of the FRS: 

- (18) 



For identified fragments, A and q are known as integer values. Therefore, 
the precision obtained on the velocity depends only on the precision on the 
magnetic rigidity. As previously discussed, it can be estimated to be in the 
order of 3.10"^, which is far better than the best time-of- flight measurements 
for heavy ions. 




(Py)i (Py)i m. 

Fig. 14. Velocity distribution for lead isotopes in the first part of tfie FRS during 
tfie 500A MeV Ph+p experiment. Only fragments ttiat kept the same charge state 
in both parts of the spectrometer are represented. The continuous lines correspond 
to the fit described in the text. 

Any error on the estimation of the charge state of fragments has a direct 
consequence on the calculated mass and therefore on the calculated momen- 
tum value. Let us consider for example a fragment {Z,A) which was assumed 
to be fully stripped in both parts of the FRS. If this fragment was in fact 



28 



an hydrogcn-like ion, its real mass is ^4 — 3 and its calculated momentum is 
related to its real one by: 



in \ (calculated) {Bp)i A — 3Z 

- {A - 3)/Z - 1 • 

= (/37)r^^^ (20) 



The velocity distributions of fragmentation residues are known to be of Gaus- 
sian shape. As it can be seen on figure 14, the misidentified fragments create 
a second component in the velocity spectra. According to the previous discus- 
sion about the fragmentation experiments, this component is very important 
for high-Z, neutron-rich fragments while it slowly vanishes for low-Z and/or 
proton-rich fragments. 

There are two ways to deal with this contribution: cither a direct estimation 
through a fit, or an estimation of the whole distribution (also by a fit) followed 
by a subtraction of the part contributed by the pollution. If the first solution 
is well suited for fragments very close to the projectile, for which the velocity 
spread is small enough for the two peaks to be distinguished, its application 
to the case of lighter fragments, for which the two peaks partially or totally 
overlap, would lead to large uncertainties. Therefore, we recommend the sec- 
ond solution, and we will now present a way to implement it in a numerically 
stable way. 

The velocity distribution of each fragment can be described as the sum of 
two Gaussian functions, the second one accounting for the pollution related 
to misidentified fragments. In order to optimize its stability, a single fit is 
applied, for each set of data associated to a given combination of Z and Aq. 
The coefficients describing the center (noted (/^0')o) and the width (noted a) 
of each Gaussian arc related to one another by simple, physics-based relations 
derived from the Morrissey systematics [17]: 

if3j)PiA)^aoiZ)+a,iZ).A (21) 



a^^^ (A) = .JPo{Z)+Pi{Z).A + P2{Z).A^ (22) 
Defining an elementary velocity distribution as: 

G{A, Z){(5-f) = A(A, Z) . e ^^^^W (23) 
the velocity function can be simply written as: 
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Fig. 15. Comparison between ratios of charge-state probabilities P(0,0)/P(1,1) and 
P(0,l)/P(l,2) as obtained from the fit of the experimental momentum distributions 
with the values calculated using the code GLOBAL (see text for details). 

f(^^^)^P^) = GiA,Z)iP^)+piZ).G{Z,A~3)i^.^.p^) (24) 



The p{Z) coefficient stands for the ratio between the probabihty of the con- 
sidered charge state and the probabihty of the charge state responsible for 
contamination. It is a free parameter of the fit. The value obtained for each 
Z can be compared to values calculated with a code dedicated to ionic charge 
states like GLOBAL [5]. As can be seen on figure 15, in the 500A MeV Pb+p 
experiment, the overall agreement was satisfactory. 

The experimental production rate I for each isotope, including the contribu- 
tion from misidentified fragments, is obtained by integration of the / function 
defined in equation 24. 

6.2 Correction of production rates for misidentified fragments 

The contribution from misidentified fragments can be subtracted iteratively: 
P(gi, q2).T{A, Z) = I{A, Z, q,, q^) - P{q, + 1, + 1).T{A - 3, Z) (25) 



As the production cross section of the lightest isotopes decreases exponentially, 
iteration can be started considering that the net production rate T is for the 
hghtest isotopes. To calculate the charge-state probability P{qi + I,q2 + 1) we 
used the GLOBAL code. We considered a 10% uncertainty for this calculation. 
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Fig. 16. Production rates for Hg isotopes obtained for different charge states combi- 
nations (left); same data, each distribution being normalized by its respective ionic 
charge state probability (right). 

This leads to a high uncertainty for the few neutron-rich isotopes, for which 
the correction reaches several tens of percent. 



6.3 Normalization of production rates 



At this step, production rates have been obtained for all the most likely charge- 
state combinations. In order to reduce the statistical uncertainty, it is appro- 
priate to consider the production rates obtained for all combinations. This 
differs from high-energy experiments in which it is generally convenient to 
disregard the non-fuUy stripped fragments. 

The simplest way to normalize all the different production rates consists in 
using calculated charge-state probabilities, by example using GLOBAL. But it 
is also possible to fit those probabilities by requiring agreement on the differ- 
ent normalized distributions. This process is illustrated in figure 16 for the Hg 
isotopes measured in the 500A MeV experiment. The agreement between the 
various distributions is very satisfactory. The use of the multiple charge-state 
combinations has two interests: a reduction of the uncertainty on the produc- 
tion rate and access to an experimental value for the charge-state probabilities. 



Further check of this method can be achieved by comparing the values obtained 
for the charge-state probabilities to GLOBAL calculations. This comparison is 
presented in figure 17 for the 500A MeV experiment. The agreement is better 
than 10% for most combinations. Nevertheless a spectacular deficiency of the 
code is observed for Li-like ions for which production is underestimated by 
GLOBAL by a factor 2.5. 

From this point, the goal of getting the cross sections can be achieved in the 
same way as in any experiment: one needs to divide the production rates by 
the number of atoms in the target and by the number of ions from the beam, 
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Fig. 17. Comparison of probabilities of charge-state combinations, plotted as the 
ratio of experimental values (see text) over GLOBAL calculations. 

and correct for the efficiency of the detectors. In the case of fragmentation 
reactions, multiple reactions play an important role in the production of light 
isotopes. This is especially true at low energy as the production cross sec- 
tions fall sharply [18]. Therefore, the contribution of these multiple reactions 
must be carefully subtracted. A dedicated method has been developed. This is 
largely self-consistent and uses a limited input from the EPAX parametriza- 
tion [19]. This method will be presented in a forthcoming paper [18]. 



7 Conclusions 



It has been demonstrated that a clean Z identification can be obtained for 
reaction products even when they are not completely stripped by increasing 
the quantity of gas of the ionization chambers with respect to the quantities 
used for high-energy experiments. 

Conclusions are not so clear regarding mass identification, which requires the 
charge-state measurement using a thick degrader. According to our calcu- 
lations, this identification is possible at 500A MeV, although the expected 
resolution is extremely close to the acceptable limit. We have shown that dur- 
ing the 500A MeV Pb-l-p experiment the obtained resolution was far worse 
than the calculated value. Using additional simulations we propose that the 
main reason for this failure is the presence of surface inhomogeneities on the 
degrader. Therefore, it may be possible that the technique developed for this 
experiment will be successfully applied in the future, when an improved de- 
grader will be available. 
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We have demonstrated that it is possible to obtain production rates with a 
high accuracy, even in the case of a failure of the charge-state measurement. 
This can be achieved through a systematic evaluation of the contribution of 
the different charge states and the suppression of the contribution of polluting 
ones. Calculations by an atomic-physics code (GLOBAL) are needed in order 
to estimate some corrections, but most of the atomic-physics data can be 
obtained in a self-consistent way from the data of the experiment. In the case 
of the 500A MeV experiment, comparison between GLOBAL calculations and 
obtained atomic-physics data confirmed the overall accuracy of this code at 
this energy, with the exception of the estimation of Li-like ions. 

If one considers the low-energy limit defined in the early papers dedicated to 
the FRS, the overall success of the 500A MeV experiment is clearly a push 
beyond this limit. Nevertheless, further decrease in energy for experiments in 
the system Pb-l-p would create further problems, as further multiphcation of 
the charge states would lead to additional complexity in the analysis. 



A Relation between straggling in energy and magnetic-rigidity 
width 

An energy loss can be written as a normalized value: 

AE = A(7i-72) (A.l) 



Inserting the unit-less version of equation 2, the previous equation writes: 

ii + {{Bp),.^y - Ji + {{Bph.^y (A.2) 



In section 4.3.1 we have shown that the contribution of (-Bp)i to the total 
uncertainty is negligible. This leads to the dispersion in {Bp)2 and the energy 
loss: 



Dividing by the energy loss and replacing the Bp.q/A terms by the Lorentz 
coefficients, one gets: 

5{^E) _ (/37), 5{Bp), 
AE ^'■71-72' (Bp), ^ ■ ' 
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One can eventually introduce the difference of magnetic rigidities in the above 
equation: 

AE y"'- 71-72 )■ (Bp), - {Bph ■ (Bp), ^ ' ' 

In the energy range considered in this paper (a few hundreds of MeV to one 
GeV per nucleon) and with the degrader thickness selected as half of the range 
of the beam, the kinematic factor (left part of the expression in equation A. 5) 
is always close to 1 (for example, 0.92 at 500A MeV or 0.95 at lA GeV). This 
could be expected, as this term may be approximated here hj P2/ < P >i^2- 
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